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Abstract. 

Both regular (the zero-momentum ghost dressing function not diverging), also named decoupling, and critical (diverging), 
also named scaling, Yang-Mills propagators solutions can be obtained by analyzing the low-momentum behaviour of the 
ghost propagator Dyson-Schwinger equation (DSE) in Landau gauge. The asymptotic expression obtained for the regular 
or decoupling ghost dressing function up to the order ff(q^) fits pretty well the low-momentum ghost propagator obtained 
through the numerical integration of the coupled gluon and ghost DSE in the PT-BFM scheme. Furthermore, when the size of 
the coupling renormalized at some scale approaches some critical value, the PT-BFM results seems to tend to the the scaling 
solution as a limiting case. This critical value of the coupling is compared with the lattice estimate for the Yang-Mills QCD 
coupling and the latter is shown to lie much above the former. 
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INTRODUCTION 

The low-momentum behaviour of the Yang-Mills propa- 
gators derived either from the tower of Dyson-Schwinger 
equations (DSE) or from Lattice simulations in Landau 
gauge has been a very interesting and hot topic for the 
last few years. It seems by now well established that, 
if we assume in the vanishing momentum limit a ghost 
dressing function behaving as F{q'^) ^ [q^)"'' and a 
gluon propagator as A((7^) ^ {q^)"^-^ (or, by following 
a notation commonly used, a gluon dressing function as 
G{q^) = q^A{q^) ^ (^^)"°), two classes of solutions may 
emerge (see, for instance, the discussion of refs. [1, 2]) 
from the DSE: (i) those, dubbed "decoupling", where 
Up = and the suppression of the ghost contribution 
to the gluon propagator DSE results in a massive gluon 
propagator (see [3, 4] and references therein); and (ii) 
those, dubbed "scaling", where ap ^ and the low- 
momentum behaviour of both gluon and ghost propaga- 
tors are related by the coupled system of DSE through 
the condition lap + OCg = implying that F^ {q^)G{q^) 
goes to a non-vanishing constant when q^ ^0 (see [5, 6] 
and references therein). 

Lattice QCD results appear to support only the mas- 
sive gluon (ac = 1) or scaling solutions (see [11] and 
references therein), and also pinching technique results 
(see, for instance, [12, 18] and references therein), re- 
fined Gribov-Zwanziger ' formalism (see [13]) or other 



' In addition, K-I. Kondo triggered very recently an interesting dis- 
cussion about the Gribov horizon condition and its implications on the 
Landau-gauge Yang-Mills infrared solutions [7. 8, 9, 10]. 



approaches like the infrared mapping of X^"^ and Yang- 
Mills theories in ref. [14] or the massive extension of the 
Fadeev-Popov action in ref. [15] appear to point to. 

In the present contribution, we will briefly review the 
work of ref. [17], which extended the previous studies 
of refs. [1, 2, 16], by the analysis of the results obtained 
by solving the coupled system of Landau gauge ghost 
and gluon propagators DSE within the framework of the 
pinching technique in the background field method [18] 
(PT-BFM) 

THE TWO KINDS OF SOLUTIONS OF 
THE GHOST PROPAGATOR 
DYSON-SCHWINGER EQUATION 

As was explained in detail in refs. [2, 16, 17], the 
low-momentum behavior for the Landau gauge ghost 
dressing function can be inferred from the analysis of 
the Dyson-Schwinger equation for the ghost propagator 
(GPDSE). That analysis is performed on a very general 
ground: one applies the MOM renormalization prescrip- 
tion, Fr{ii^) = ^^Ar{h^) = 1, where ji^ is the subtrac- 
tion point, chooses for the ghost-gluon vertex, 

rf'{-q,k;q-k) = igof"''- {qvHi{q,k) 

+ {q~k)yH2iq,k)) (1) 

to apply this MOM prescription in Taylor kinemat- 
ics {i.e. with a vanishing incoming ghost momentum) 
and assumes the non-renormalizable bare ghost-gluon 
form factor, H\{q,k) — H\, to be constant in the low- 
momentum regime for the incoming ghost. Then, the low 



momentum-behaviour of the ghost dressing function is 
supposed to be well described by 
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and that of the gluon propagator by 
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and this, after solving asymptotically the GPDSE, finally 
left us with: 
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It should be understood that the subtraction momentum 
for all the renormalization quantities is ji^. The case 
ap j^O leads to the so-called scahng solution, where the 
low-momentum behavior of the massive gluon propaga- 
tor forces the ghost dressing function to diverge at low- 
momentum through the scahng condition: 2a^ + ofe = 
(ocg = 1 is the power exponent when dealing with a mas- 
sive gluon propagator). As this scaling condition is veri- 
fied, the perturbative strong coupling defined in this Tay- 
lor scheme [19], Ut = g^/{4n), has to reach a constant 
at zero-momentum. 
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as can be obtained from Eqs.(3,5). The case ttf = cor- 
responds to the so-called decoupling solution, where the 
zero-momentum ghost dressing function reaches a non- 
zero finite value and eq. (5) provides us with the first 
asymptotic corrections to this leading constant. This sub- 
leading correction is controlled by the zero-momentum 
value of the couphng defined in eq. (6), which is an ex- 
tension of the non-perturbative effective charge defini- 
tion from the gluon propagator [20] to the Taylor ghost- 
gluon coupling [21]. 



TABLE 1. Gluon masses and the zero- 
momentum non-perturbative effective 
charges, taken from ref. [17] and ob- 
tained as discussed in the text. 
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0.15 


0.24 


0.37 


0.16 


0.30 


0.39 


0.17 


0.41 


0.43 



NUIVIERICAL RESULTS FROIVI 
COUPLED PT-BFJVI DSE'S 

In ref. [17], the solutions of the coupled DSE system in 
the PT-BFM scheme (with H\ = I for the ghost-gluon 
vertex), numerically integrated for many values of the 
coupling at the renormalization point /x^ as a boundary 
condition, were studied in the hght of the analytical 
results above presented. Here we will shortly discuss the 
results of this work. 



The "regular" or "decoupling" solutions 

The numerical results of the PT-BFM coupled DSEs 
were shown in ref. [17] to behave asymptotically as 
eq. (4) predicts for the decouphng DSE solutions. In- 
deed, as the gluon propagator solutions in the PT- 
BFM scheme result to behave as massive ones, the 
eqs. (3,5) must account for the low-momentum be- 
haviour of both gluon propagator and ghost dressing 
function with H\ = I and aj (0) given by eq. (6), with 
ar(jLl^) = gji{jJ.^) I (An) being fixed, as a boundary con- 
dition for the numerical integration of the coupled DSE 
for each particular solution of the family (see tab. 1). Fur- 
thermore, the zero-momentum values of the ghost dress- 
ing function, /^^(O) and of the gluon propagator, Aff(O), 
can be taken from the numerical solutions of the DSE 
(for any value of the a{\l = lOGeV)). These altoghether 
with the gluon masses obtained by the fit of eq. (3) to 
the numerical DSE gluon propatator solutions (see the 
left plot in fig. 1, for a(ju) = 0.16, and the results for 
a(ju) = 0.15,0.16,0.17 in tab. 1, taken from ref. [17]), 
provide us with all the ingredients to evaluate, with no 
unknown parameter, eq. (5). 

Indeed, the expression given by eq. (5) can be succes- 
fully applied to describe the solutions all over the range 
of couphng values, a(/i), at = 10 GeV (provided that 
they are not very close of the critical coupling that will 
be defined in the next subsection). This can be seen, for 
instance, for a = 0.16,0.18, in the right plots of fig. 1 
and it was also shown for a = 0. 15 in ref. [17]. 
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FIGURE 1. Gluon propagator (left) and ghost dressing function (riglit) after tiie numerical integration of the coupled DSE 
system for a{fl = lOGeV) = 0.16,0.17 taken from [17] . The curves for the best fits to gluon propagator and ghost propagator 
data explained in the text appear as red dotted lines, the same for the black dotted line in the lefthand plot but retaining only the 
logarithmic leading term for the asymptotic ghost dressing function by dropping the —11/6 away. 



The "critical" limit 

There also appeared to be a critical value of the cou- 
pling, acrit = o:(/x^) ~ 0.182 with — lOGev, above 
which the coupled DSE system does not converge any 
longer to a solution. As a matter of fact, we know from 
refs.[2, 17] that the scaling solution implies for the cou- 
pling 
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where B{^^) and A(/i^) defined by Eqs. (3,4). This is 
also shown in ref. [1], where only the ghost propagator 
DSE with the kernel for the gluon loop integral is ob- 
tained from gluon propagator lattice data. In the analy- 



sis of ref. [1], a ghost dressing function solution diverg- 
ing at vanishing momentum appears to exist and verifies 
eqs. (4,8), while regular or decoupling solutions exist for 
any a < acrit- In ref. [17], a more complete analysis is 
performed: first by studying the solutions for many dif- 
ferent values of the coupling, a ~ a{ji^), of a coupled 
DSE system; and then by showing that the ghost dressing 
function at vanishing momentum, F{0,jj.^), is described 
by the following power behaviour. 



^(0) - (acrit -a(Ai')) 



(9) 



where is a critical exponent (depending pre- 

summably on the renormalization point, ji^), supposed 
to be positive and to govern the transition from decou- 
pling (a < acrit) to the scaling (a = acrit) solutions; 



and where we let Ocrit be a free parameter to be fit- 
ted by requiring the best Hnear correlation for log[F(0)] 
in terms of log[acrit — cc]. In doing so, the best corre- 
lation coefficient is 0.9997 for Ofcrit = 0.1822, which is 
pretty close to the critical value of the coupling above 
which the coupled DSE system does not converge any 
more, and k{ij,^) = 0.0854(6). This can be seen in fig. 2, 
where the log-log plot of Fr{0) in terms of a^nt ~ CC 
is shown and the linear behaviour with negative slope 
corresponding to the best correlation coefficient strik- 
ingly indicates a zero-momentum ghost propagator di- 
verging as a ^ acrit- Nevertheless, no critical or scal- 
ing solution appears for the coupled DSE system in the 
PT-BFM, although the decoupling solutions obtained for 
any a < Ofcrit = 0.1822 seem to approach the behaviour 
of a scaling one when a Ucm- This is again well un- 
derstood in ref. [17], where the gluon propagators ob- 
tained from the coupled DSE system in PT-BFM were 
also found to obey the same critical behaviour pattem 
as the ghost propagator, when approaching the critical 
value of the coupUng. Indeed, when approaching the crit- 
ical value of the coupling, the gluon propagators obtained 
from the coupled DSE system in PT-BFM must be also 
thought to obey the same critical behaviour pattem as 
the ghost propagator. In the PT-BFM, the value at zero- 
momentum being fixed by construction [4], one should 
expect that, instead of decreasing, the gluon propagator 
obtained for coupUngs near to the critical value increases 
for low momenta: the more one approaches the critical 
coupling the more it has to increase. This is indeed the 
case, as can be seen in fig. 2(b). This also implies that, 
near the critical value, the low momentum propagator 
does not obey eq. (3) and that consequently eq. (5) does 
not work any longer to describe the low momentum ghost 
propagator. 



of flavours. One can use instead the available perturba- 
tive four-loop formula describing the running of Ihe cou- 
pling in Taylor-scheme to estimate Aqcd in this particu- 
lar scheme, then perform the conversion to MS (see for 
instance eqs. (22,23) of the first reference in [19]) and 
thus obtain the value quoted in tab. 2. Of course, it would 
be again meaningless to compare this last value with the 
one for A-^ that can be obtained from the PDG value 
for aj35(Mz), also quoted in tab. 2, but we can refer the 
comparison to the lattice Yang-Mills determinations of 
the same parameter, as for instance the two of them in- 
cluded in tab. 2. It should be noted that the procedures 
for the lattice determination of mainly work in the 
UV domain, where IR sources of uncertainties as the Gri- 
bov ambiguity or volume effects are indeed negligible ^. 
Thus, the lattice estimates of A^ appear to lie clearly 
below this critical hmit for the PT-BFM DSE in pure 
Yang-Mills. This last results appears to indicate that the 
critical solution is not the one chosen by the zero-flavour 
world. However, as no quark flavour loops effect have 
been incorporated in our DSE analysis, we cannot yet 
neither compare with the physical strong coupling nor 
conclude whether the critical limit can be allowed in the 
"real world". 

TABLE 2. The critical value of Aj^ in pure Yang-Mills 
inferred from acrit = 0.1822 (first column), lattice estimates 
for Yang-Mills taken from literature (second and third 
columns) and the one obtained from the PDG value of 
O^(Mz) by applying a four-loop perturbative formula for the 
running of ct^ with Nj — 5. 



MS.crit 








434 MeV | 238(19) Me V | 244(8) MeV | 213(9) MeV 



CONCLUSIONS 



The critical value and the QCD coupling in 
the "real world" 

Finally, one can pay attention to the critical value of 
the couphng, acrit = 0.1822, and try to make a compar- 
ison with the physical strong coupling values in order 
get some idea of whether the current data can exclude 
or not this critical behaviour. Although the experimen- 
tal PDG world average of the strong coupling in the 
MS scheme, c^(Mz) = 0.1184(7) [22], can be prop- 
agated from the Z" boson mass down to ;U = 10 GeV 
to give aMs(10 GeV) = 0.179(2), that incidentally lies 
on the right ballpark of the above critical value, such a 
comparison is meaningless because our coupling corre- 
sponds to one in MOM Taylor- scheme for zero number 



The ghost propagator DSE, with the only assumption of 
taking Hi{q,k) from the ghost-gluon vertex in eq. (1) to 
be constant in the infrared domain of q, can be exploited 
to look into the low-momentum behaviour of the ghost 
propagator. The two classes of solutions named "decou- 
phng" and "scaling" can be indentified and shown to de- 
pend on whether the ghost dressing function achieves a 
finite non-zero constant (ay =0) at vanishing momen- 
tum or not (ap ^ 0). The solutions appear to be dialed by 
the size of the couphng at the renormalization momen- 
tum which plays the role of a boundary condition for the 



^ As a matter of fact, there are unquenched lattice determinations with 
Nf = 5 staggered fermions for the strong coupling [23] which are pretty 
consistent with the PDG value. This can be taken as a good indication 
in favour of the robustness of the lattice determinations of A-nw. 



A(q",H") 




FIGURE 2. (a) Log-log plot of the zero-momentum values of the ghost dressing function, obtained by the numerical integration 
of the coupled DSE system in the PT-BFM scheme, in terms of oCaU — a. a — a{ll = lOGeV), the value of the coupling at the 
renormalization momentum, is an initial condition for the integration; while 0!ei-it is fixed to be 0.1822, as explained in the text, by 
requiring the best linear correlation, (b) Gluon propagator solutions in terms of for the same coupled DSE system for different 
values of a(/x — lOGeV), all very close to the critical value, ranging from 0.18 to 0.1817 . 



DSE integration. The low-momentum behaviour of the 
decoupUng solutions results to be regulated by the gluon 
propagator mass and by a regularization-independent di- 
mensionless quantity that appears to be the effective 
charge defined from the Taylor-scheme ghost-gluon ver- 
tex at zero momentum. 

In this note, we have shortly discussed the results 
of ref. [17] where the solutions of coupled ghost and 
gluon propagator DSE in the PT-BFM scheme were 
studied and demonstrated that the asymptotic decou- 
pling formula (ap = 0) successfully describes the low- 
momentum ghost propagator The model applied for the 
massive gluon propagator is also verified to give prop- 
erly account of the gluon solution, at least for momenta 
below 1 GeV (and for a coupling not very close to the 
critical point). Although we argued that a massive gluon 
propagator implies that the ghost dressing function takes 
a non-zero finite value at vanishing momentum, we also 
show that the zero-momentum ghost dressing function 
tends to diverge when the value of the coupling dial- 
ing the solutions approaches some critical value. Such 
a divergent behaviour at the critical coupling seems to 
be the expected one for a scaling solution (where, if 
the gluon is massive, ap — —1/2). If we consider the 
zero-momentum value of the ghost dressing function 
as some sort of "order parameter" indicating whether 
the ghost propagator low-momentum behaviour is sup- 
pressed (ttf = and finite ghost dressing function) or it 
is enhanced (a^ < and divergent ghost dressing func- 
tion), the strength of the coupling computed at some 
renormalization point seems to control some sort of tran- 



sition from the suppressed to the en/zancec/ phases for the 
ghost propagator DSE solutions in the PT-BFM scheme. 
The last only takes place as some critical value of the 
coupling is reached. Neverteless, it can be proven that, 
as far as the gluon is massive, the scaling behaviour for 
the Yang-Mills propagators appear not to be a solution 
but an unattainable limiting case for the PT-BFM DSE 
solutions. 

Finally, the critical value for the coupling is shown to 
lie clearly much above the estimate of the QCD coupling 
in pure Yang-Mills computed from lattice QCD. This of 
course agrees with the fact that the current large-volume 
quenched lattice results for ghost and gluon Green func- 
tions clearly behave as expected for the decoupling solu- 
tions. 
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